Collective excitations of Bose-Einstein condensed gases at finite temperatures 
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We have applied the Popov version of the Hartree-Fock- 
Bogoliubov (HFB) approximation to calculate the finite- 
temperature excitation spectrum of a Bose-Einstein conden- 
sate (BEC) of 87 Rb atoms. For lower values of the tempera- 
ture, we find excellent agreement with recently-published ex- 
perimental data for the JILA TOP trap. In contrast to recent 
comparison of the results of HFB-Popov theory with exper- 
imental condensate fractions and specific heats, there is dis- 
agreement of the theoretical and recent experimental results 
near the BEC phase transition temperature. 

PACS Numbers: 3.75.FL 67.40.Db, 67.90. +Z 



Laboratory realizations of gaseous Bose-Einstein con- 
densates (BECs) have prompted vigorous experi- 
mental Q investigations of the temperature-dependent 
properties of these mesoscopic quantum systems. The 
current theoretical interest in such condensates derives, 
in part, from the fact that experimental tests of many- 
body theories that are thought to apply to BEC — a 
phenomenon occuring in many areas of physics — can be 
performed for the first time. An accurate theory of such 
systems is therefore of fundamental interest, and will also 
have practical applications. In this Letter, we explore the 
limits of validity of the simplest temperature-dependent 
mean-field theory — a simplified version of the Hartree- 
Fock-Bogoliubov (HFB) approximation originally intro- 
duced by Popov |g[ — by presenting the first comparison 
with experiment of this theory's predictions for conden- 
sate excitation spectra at T > 0. 

Condensate properties predicted by such theories in- 
clude condensate and thermal-atom spatial density pro- 
files, condensate fractions, specific heats, and excitation 
frequencies. Of these properties, excitation spectra pro- 
vide the most sensitive test of the applicability of com- 
peting theories since the other quantities listed depend 
on sums over states and are thus insensitive to small er- 
rors in the excitation spectrum. One example of this 
can be found in the approach of zero-temperature exci- 



tation frequencies to the Thomas-Fermi limit (i.e., the 
limit No — » oo, where Nq is the number of condensate 
atoms) as Nq increases. For 87 Rb condensates confined 
in the JILA TOP (time-averaged orbiting potential) trap 
where No > 4000 atoms, the Thomas-Fermi predictions 
for spatial density profiles are extremely accurate while 
low-order excitation frequencies can differ by as much as 
10%. 

There is good a priori reason to expect that the 
HFB-Popov theory should provide good predictions of 
experimental T-dependent collective excitation frequen- 
cies. Measurements of zero-temperature frequencies || 
exhibited excellent agreement with the predictions of 
zero-temperature, mean field theory ffl]. Furthermore, 
semiclassical variants of the HFB-Popov theory have 
exhibited excellent agreement with experiment for In- 
dependent condensate fractions and specific heats for 
temperatures up to near T c ||. The HFB-Popov the- 
ory is a finite-temperature extension of mean-field the- 
ory which provides self-consistent treatment of the con- 
densed and thermal components of the gas and which 
should describe the linear response of the condensate to 
small-amplitude mechanical disturbances ||. 

Although the HFB-Popov equations have been derived 
elsewhere HfLO), we shall briefly state the physics be- 
hind the basic equations here. The confined Bose gas is 
portrayed as a thermodynamic equilibrium system un- 
der the grand canonical ensemble whose thermodynamic 
variables are N, the total number of trapped atoms, T, 
the absolute temperature, and either No or [i, the chem- 
ical potential. The system hamiltonian has the form 

K = H-fiN = I dnp\r)(H - n)i>(r) 



where 4>(y) is the Bose field operator which annhililates 
an atom at position r and Hq = "jff V 2 + Vtrap( r ) is 
the bare trap Hamiltonian. For the system treated here, 
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the trap potential Vtra P (r) = M {to 2 p 2 + <^ 2 z 2 ) /2, where 
M is the atomic mass, and lo p and lo z are the radial 
and axial trap frequencies, respectively. The quantity 
Uo = 4:Trh 2 a/M is a measure of the interaction strength 
between atoms, with a being the scattering length for 
zero-energy binary atomic collisions, taken to be 109ao 
for 87 Rb where ao is the Bohr radius. 

The Bosc field operator is decomposed into a c-number 
condensate wave function plus an operator describing the 

A 1/2 ~ 

non-condensate part: -0( r ) = N ' "M 1 ") + V'( r ) an d in- 
serted into Eq. (|l|). When terms cubic and quartic in 
ip(r) are treated within the mean-field approximation 
the grand canonical hamiltonian reduces to a sum of four 
terms: K = Kq+Ki+k\+K2- The first term, Kq is a c- 
number, the second and third terms are linear in ip(r) and 
ip (y) and the last term is quadratic in these quantities. 
It is easy to show that the linear terms vanishes iden- 
tically if <fi(r) satisfies the generalized Gross-Pitaevskii 
(GP) equation 

{h + U Q [N |0(r)| 2 + 2n(r)] } 0(r) - t*P(r), (2) 

where n(r) is the density of non-condensate atoms. Note 
that the condensate wavefunction is normalized to unity. 
The term K2 has the form 

K 2 = [ drft (r)&p(r) + / dr (0) 2 (r)^ (r)^ (r) 



where Nj — (eP E: > — l) 1 is the Bose-Einstein factor, 
and (3 = (fcs? 1 ) 1 with /cb the Boltzmann constant. The 
total number of trapped atoms, TV, is given by 



+ AW, I rfr(r) 2 (r) ^ (r) ^ (r) (3) 

where C = H + 2Uon(r) — \i and n(r) = iVo |0(r)| 2 + n(r) 
is the total trapped-atom density. The term K2 can be 
diagonalized by the following Bogoliubov transformation 

^(r) = J2( u j( v ) a j + v *j » 

3 

^t( r ) = ^( u *( r ) a t +Ui ( r ) a .) . (4) 
3 

if the quasi-particle amplitudes Uj(r) and Vj(r), satisfy 
the coupled HFB-Popov equations: 

£ Uj {r) + N U |^(r)| 2 Vj (r) = E jUj (r) 

£ Vj (r) + N U Q |0(r) | 2 Uj (r) = (r) , (5) 

The aj and a| are quasi-particle annihilation and cre- 
ation operators that satisfy the usual Bose commutation 
relations. 

The density of the thermal component of the gas 
n(r) = (ip'ip) and thus can be written in terms of the 
quasi-particle amplitudes as 

n(r) = J2 { [KMI 2 + \ v ^)f} N 3 + - (6) 
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dr n(r) = Nq + dr h(r) 



(7) 



Our version of Eq. (§) differs fr om that of Rcf. H via a 
sign change in the definition of Vj(r). 

Equations (|^), (||), (^|), and (^) form a closed system of 
equations that we have referred to as the "HFB-Popov" 
equations. We have numerically solved these equations 
under conditions appropriate to 87 Rb atoms confined in 
the JILA TOP trap jl3|]. We choose our state variables 
to be {T, /i, N}, fix T and /x, and then determine N by 
solving the HFB-Popov equations. This is equivalent to 
the alternative triple of state variables {T, No, N}, since 
there is a one-to-one relationship between A^o and /i. 

We have solved the HFB-Popov equations by an iter- 
ative procedure, each cycle of the iteration consisting of 
two steps. In the first step of each cycle, we solve Eq. 
(g) for new values of ^(r) and No with a basis-set ap- 
proach as described previously [jlif using n(r) obtained 
in the previous cycle. In the second step, we solve Eqs. 
(||) using n(r) from the previous cycle, and the newly 
generated values of 4>{r) and TVg. With the quasi-particle 
amplitudes expanded in the trap basis, Eqs. (||) yield a 
generalized matrix eigenvalue problem for the basis-set 
coefficients. We recast the generalized matrix eigenvalue 
problem by using a decoupling transformation consist- 
ing of taking the sum and difference of Eqs. (||). This 
transformation is equivalent to that of Hutchinson et 
al. , except that it is expressed in terms of basis set 
expansion coefficients. Completion of this step yields the 
{uj(r), Vj(r)}, and Ej which are used in Eq. (||) to up- 
date n(r). Equation (^) then updates the total number 
of trapped atoms, N. Convergence is reached when the 
change in TV from one cycle to the next is smaller than a 
specified tolerance. To obtain converged results at high- 
temperatures, we add a correction to the total number 
of atoms, N, at each iteration cycle. High-energy quasi- 
particle eigenfunctions have negligible overlap with the 
condensate wave function, so their presence in the ther- 
mal sum of Eq.(||) does not significantly modify the low- 
lying excitation frequencies, but does contribute to the 
value of N. 

We have checked the accuracy of our numerical work 
by writing two independent codes, which produce iden- 
tical answers. The ideal gas result is recovered when we 
set a — 0, and we have reproduced the results reported 
in Ref. |Q. We now discuss the comparison of this ap- 
proach with experiment. 

Figure [TJ compares the experimental Jl3j excitation 
spectrum of 87 Rb in the JILA TOP trap, vs. our HFB- 
Popov results for the m = and m = 2 modes. The ab- 
scissa is the scaled temperature T' = T/T (N,lu) where 
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T (N,uj) = hu/k B [7V/C(3)] 1/3 is the theoretical tran- 
sition temperature for an ideal, trapped Bose gas and 

1/3 

uj = (ujpUj z j ■ The ordinate is the excitation frequency 
expressed in units of u> p . Our results were obtained us- 
ing the experimental value of T, and a value of /i that 
yielded the experimentally determined value of N. Thus, 
as for our previous treatment of zero-temperature exci- 
tation spectra Q, this calculation contains no adjustable 
parameters. The agreement between theory and experi- 
ment is very good (on the order of 5%) for low and in- 
termediate temperatures (T" < 0.65). It should be noted 
that, as depicted in Fig. ||, the high end of this tem- 
perature range corresponds to a non-condensate fraction 
of about 50%. However, as the temperature increases, 
the HFB-Popov excitation frequencies diverge from the 
experimental data. This feature of the comparison holds 
true for both m = and m = 2 modes. 

The behavior of the calculated excitation frequencies 
can be understood in a simple way. The HFB-Popov 
equations determine the equation of state for the state 
variables {N,Nq,T}, so that, given the values of N and 
T, a unique Nq is determined. For fixed N, this relation- 
ship generates the condensate fraction Nq/N as a func- 
tion of T, which is shown in Fig. || for the JILA TOP trap 
with N — 2000. One can easily predict the temperature- 
dependent mode frequencies for the N — 2000 system 
by finding the number of condensate atoms, No, from 
Fig. ||, and then determining the zero-temperature exci- 
tation frequency of a condensate with No atoms, which is 
a much simpler calculation. Figure [| shows a comparison 
of the three lowest frequencies numerically determined 
from the HFB-Popov equations, and by the equivalent 
T = method just described. The agreement of these 
two approaches is very good over nearly the entire tem- 
perature range. The two solid curves in Fig. |]are the fre- 
quencies determined by the same procedure except that 
the number of condensate atoms was taken from exper- 
iment. In short, the principal effect of finite tempera- 
ture on the HFB-Popov excitation spectra is largely an 
effect of condensate depletion: the dynamics of the finite- 
temperature condensate are essentially the same as those 
of a zero-temperature condensate with the same value of 
Nq. This is consistent with earlier calculations [ p||l6| 
of the speed of sound in a homogeneous Bose-condensed 
gas, which found that its temperature dependence was 
effectively a condensate-density dependence. We discuss 
this result In a broader context in a separate paper [jl7| , 
in which we show that HFB-Popov results can be re- 
produced quantititatively by a much simpler "two-gas" 
model: the condensate gas, which is described by the 
zero-temperature GP equation; and the thermal cloud, 
which is described as an ideal Bose gas in an effective 
potential created by the condensate. This effective po- 
tential repels the thermal gas from the condensate gas, 
which results in the essential independence of tempera- 



ture of all condensate properties except Nq. Application 
of simple quantum statistical mechanics to this model 
can generate the full phase diagram of Fig. ^ directly. 

As Fig. [j] clearly shows, the HFB-Popov solutions re- 
produce the experimental results quite well when T < 
0.65 To, but fail at higher temperatures. The HFB-Popov 
formalism is biased toward description of the condensate, 
as it represents the condensate excitations as taking place 
in a static thermal cloud. This results in at least one 
minor failure of the approach, which is weakly visible 
in Fig. |^ as a deviation of the m = 1 mode frequency 
from unity near To, in violation of the generalized Kohn 
theorem for parabolic confinement ]l2,|l£|. This mode 
should correspond to a rigid oscillation of the complete 
iV-atom system, and the deviation of its frequency from 
unity results from the HFB-Popov approximation hold- 
ing the thermal component fixed and allowing only the 
condensate to oscillate. In experiments of the type dis- 
cussed here, however, the thermal component and con- 
densate must both be driven by the modulation of their 
common confining potential. For other, non-rigid, oscil- 
lations we may thus also expect thermal and condensate 
modes to be coupled in general. Thus, HFB-Popov fre- 
quencies will only correspond to the experimental val- 
ues if the condensate response to mechanical disturbance 
does not induce modulations of the thermal density, so 
that there is no back-action of the thermal cloud on the 
condensate motion. A more general theory that accounts 
for such condensate-cloud interactions has recently been 
outlined [JOt, but remains to be implemented. The self- 
consistent inclusion of pair terms, which are negelected in 
the Popov approximation, may also be important for cap- 
turing multiple-collision effects in the theory of trapped 
BECs p(| . The effect of such terms can be very marked 
close to the transition temperature El]]. 

In conclusion, we have delineated the region of validity 
of the Popov version of finite temperature HFB theory 
by comparing it directly with the results of recent exper- 
iments. Good agreement is obtained for condensate frac- 
tions from unity down to about 0.5, so the HFB-Popov is 
apparently correctly describing finite-temperature phe- 
nomena in a nontrivial regime. This comparison confirms 
the critical role of evaporatively cooled gases in estab- 
lishing proper finite-temperature field theories of Bose- 
Einstein condensation, and shows that there is still work 
needed to establish satisfactory agreement between the- 
ory and experiment for cases of small condensate fraction. 
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FIG. 1. The experimental, temperature dependent excita- 
tion spectrum in the JILA TOP trap (filled circles) versus the 
HFB-Popov predictions for the m = mode (top, labeled by 
"+") and the m = 2 mode (bottom, labeled by "x"). The 
solid curves are excitation frequencies for a zero -temperature 
condensate having the same number of condensate atoms as 
the experimental condensate in the finite-T cloud. 
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FIG. 2. A plot of the condensate fraction as a function of 
T" (solid curve) for the JILA TOP trap in which N is fixed at 
2000 atoms. The same quantity is shown for the ideal gas in 
the thermodynamic limit (dotted curve) and for 2000 atoms 
(dashed curve) for comparison. 
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FIG. 3. HFB-Popov excitation frequencies (filled circles) 
for the m = (top), m — 2 (middle), and the m = 1 
modes (bottom) for a cold-atom cloud having N = 2000 
atoms. Overlaid (solid lines) are the frequencies for a 
zero-temperature system with the same number No of con- 
densate atoms as in the finite-temperature system. 
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